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ABSTRACT 

We calculate the luminosity and energy spectrum of the neutrino emission 
from electron-positron pair annihilation during the collapse of a supermassive 
star (M ^ 5 x 1O 4 M ). We then estimate the cumulative flux and energy 
spectrum of the resulting neutrino background as a function of the abundance 
and redshift of supermassive stars and the efficiency of these objects in 
converting gravitational energy into neutrino energy. We estimate the expected 
signal in some of the new generation of astrophysical neutrino detectors from 
both a cumulative background of supermassive stars and single collapse events 
associated with these objects. 

Subject headings: elementary particles: neutrinos - cosmology: observations and 
theory 
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1. Introduction 

In this paper we examine the physics of neutrino emission in the collapse of supermassive 
stars. We also comment on the prospects for future terrestrial neutrino detectors to obtain 
a signal from these objects, and we discuss the possible consequences of such a signal for 
cosmology. By "supermassive stars" we mean stars so massive that they collapse on the 
general relativistic Feynman-Chandrasekhar instability. This will imply masses for these 
objects M ^ 5 x 1O 4 M (see for example Fuller, Woosley, & Weaver 1986; thereafter 
FWW). 

Though there is no direct evidence for the existence of these objects, we note that there 
is overwhelming evidence for the existence of supermassive black holes both associated with 
quasars and Active Galactic Nuclei (AGN) at high redshift and almost every galaxy-sized 
structure examined appropriately by the Hubble Space Telescope (HST) (van der Marel et 
al. 1997). In turn, Begelman and Rees (1978) have shown that a supermassive star could 
result as an intermediate stage in the collapse of a relativistic star cluster to a black hole. 

Alternatively, supermassive stars could have formed out of primordial clouds at high 
redshifts in which cooling was not as efficient as in clouds contaminated with metals (Hoyle 
k Fowler 1963; Bond, Arnett, & Carr 1984; FWW; McLaughlin & Fuller 1996). The typical 
baryonic Jean's mass at high redshift can be ~ 10 5 M Q (Peebles & Dicke 1968; Tegmark 
et al. 1997), but we do not know whether a cloud of this size would fragment into many 
pieces and form stars of smaller masses or collapse directly to form a large object. Given 
the relatively crude understanding we currently possess on star formation, it is interesting 
to explore the observational signatures of supermassive stars, one of which could be their 
neutrino emission during core collapse. (The other telltale signs of the existence and 
evolution history of supermassive stars could be nucleosynthesis products of hot hydrogen 
burning (the rp-process, Wallace & Woosley 1981); greatly enhanced local helium and/or 
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deuterium abundances (Woosley 1977; Fuller & Shi 1997); or effects of large black holes.) 

It would be very significant for our understanding of galaxy formation and cosmology 
if there were to be a new neutrino or nucleosynthesis probe of the epoch of galaxy/quasar 
formation at redshifts z ~1 to 5. 

There are significant differences between neutrino emission in supermassive stars 
and in core collapse supernovae. Neutrinos are trapped and thermalized in ordinary 
supernovae. Because of the different depths inside the supernova core where the various 
neutrino species thermally decouple, the average neutrino energies satisfy the hierarchy 
E Vt ~ E Pt ps E Vfi ~ E P/1 > E Pe > E Ve . By contrast, in supermassive stars neutrinos are 
produced principally by annihilation of thermal e + e~-pairs. They can escape freely from 
the core and so their energy spectra are not thermal and different neutrino species therefore 
will have similar energy spectra. 

The only known limit so far on the neutrino background from supermassive stars comes 
from the consideration that probably no more than ~ (9(0.1) of all baryons could ever have 
been in these objects, otherwise there would probably be too many ~ 1O 4 M to 10 6 M© 
relic black holes. This limit is much less stringent than the limit on the supernova neutrino 
background. This supernova background limit can be obtained because the past and present 
supernova rates are subject to a very tight metallicity production constraint (Totani, Sato 
& Yoshii 1996; Hartmann & Woosley 1997; Malaney 1997). The abundance of supermassive 
stars may not be subject to similar metallicity concerns because these objects do not 
necessarily expel significant amounts of metals into the interstellar medium. However, 
while the supernova neutrino background has contributions from recent supernovae, the 
neutrino background from supermassive stars (if any) probably took its form at a higher 
redshift, and therefore suffers a fair amount of redshift in energy. Consequently, although 
the relic neutrino flux of supermassive stars could be higher than the flux of relic supernova 
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neutrinos, it has fewer neutrinos with energies in the range ^10 MeV. 

Numerical calculation of neutrino emissivity from e + e _ annihilation have been carried 
out previously (Schinder et al. 1987; Itoh et al. 1989). In this paper we apply these 
results to the calculation of the neutrino luminosity and energy spectra in the collapse 
of supermassive stars. We then proceed to calculate the flux and energy spectrum of the 
neutrino background from a putative population of supermassive stars as a function of their 
abundance and redshift, and their efficiency in converting gravitational energy into neutrino 
energy. We estimate the event rate of this neutrino background in several new-generation 
neutrino detectors, including Super Kamiokande (Super K). We will also consider the 
expected event rate of a single neutrino burst from the collapse of a single supermassive 
star. 

2. The Neutrino Luminosity from the Collapse of a Supermassive Star 

Fuller, Woosley, and Weaver (1986) have discussed the evolution and general relativistic 
instability of supermassive stars. In that work it was shown that these objects will most 
likely collapse into black holes unless the centrifugal force resulting from rapid rotation is 
strong enough to compensate the build-up of infall kinetic energy. During the collapse, only 
part of the star will plunge through the event horizon and become a "prompt" black hole. 
This is because the prodigious thermal neutrino pair emission will render the collapse of a 
nonrotating supermassive star non-homologous (FWW; see also Goldreich & Weber 1980). 
A simple calculation shows that 



where superscript "init" always refers to quantities in the initial pre-collapse configuration, 
and superscript "HC" refers to quantities of the homologous core during the collapse. Here 
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M 5 is the stellar mass in units of 10 5 M Q ; g is the statistical weight of relativistic particles, 
with g mit pa 2 and g uc ~ 2 + (7/8) x 4 = 5.5; and S is the entropy per baryon. In the initial 
nearly isentropic configuration (assuming an index n = 3 polytropic structure), we have 

where fi is the mean molecular weight (/i ~ 0.59 for "primordial composition" of 75% 
hydrogen and 25% helium by mass). In equation (0) we have assumed that g = 2 and 
that all of the entropy is contributed by photons. This is an excellent approximation for a 
supermassive star near its general relativistic instability point. 

As an example, if the entropy per baryon is reduced by a factor of 2.5 as a result of 
neutrino emission during the collapse, i.e., S nc /S imt pa 0.4, then the final homologous core 
mass will be about 10% of the initial stellar mass, Mf° /M™* pa 0.1. 

The total Newtonian gravitational binding energy of a homologous core with a mass 
M|* c is crudely ~ E s ~ 10 59 M| iC erg. If there is no strong magnetic field present, most of 
this energy will be trapped inside the black hole, radiated through gravitational waves, or 
released by the neutrino emission prior to trapped surface formation. The characteristic core 
radius near the black hole formation point is the Schwarzschild radius r s pa 3 x 10 10 M;P C cm. 
The characteristic duration of the collapse is the dynamic time t s ~ M^ c sec, or longer if 
rotation or magnetic fields hold up the collapse. 

It has been shown that the neutrino energy loss rate per unit volume, Q, as a result of 
electron-positron annihilation, has the simple form (Schinder et al. 1987; Itoh et al. 1989) 

Q pa 4 x 10 15 T 9 9 erg cm" 3 s -3 , (3) 

where Tg = T/10 9 K. This equation is valid so long as the temperature T ^ m e (the 
electron rest mass) and as long as the density p is low enough to ensure non-degeneracy 
(p ^ m p T 3 /h 3 c 3 ~ 10 8 (T/1 MeV) 3 gem -3 where m p is the proton rest mass). Both 
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conditions are satisfied during the collapse of a supermassive star. The neutrino luminosity 
from a supermassive star is then 

rR 

L u ~ ei e 2 / 4vrr 2 Q dr, (4) 
Jo 

where R is the radius of the homologous core, e\ is a factor accounting for the travel time 
difference for neutrinos coming out from different depths inside the core, and e 2 represents 
the effect of gravitational redshift. 

During a homologous collapse, the density p has a self-similar profile of polytropic 
form: 

P = pMO, (5) 

where p c is the central density of the star, and 6% is the index n = 3 Lane-Emden function 
and £ = r/a is a dimensionless length measure. The total pressure at any point in the star 
can be cast in the index n = 3 polytropic form P = Kp 4 ^. Here the pressure constant can 
be expressed as 

where Na is Avogadro's number and where g s is the relativistic particle statistical weight 
entering into the entropy per baryon S ~ (2tt 2 g s T 3 /45) / p Na- For all our considerations in 
this paper we can safely take g s = g. If we denote the density in units of 10 3 gcm~ 3 as p 3 , 
then the pressure in units of MeV 4 will be P = K 3 p^ 3 , with K 3 being 

K 3 « (6.8 x 10- MeV 4 ) (^) 1/3 ^ 3 (l + ^i) (7) 
where ^ioo is the entropy per baryon S in units of 100 Boltzman's constant. 
We can estimate similarly the dimensionless length conversion factor, 

„= IK tfV> pr V» „ ( 8 .2 x 10'° cm) (HZ!)" 6 sS (, + MLf» (-f^Y m (8) 
v 7 ^ v ^ s 1 v /i^ioo 7 v 10 3 gcm- 3 / 

where, in terms of Newton's constant G, the Planck mass is m p i = G^ 1 ^ 2 . 
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Because the entropy per baryon Si 00 is roughly constant through out the radiation- 
dominated core at any time, we can estimate that T| / p 3 = 0. 35*100 ~ constant (FWW). 
Therefore, 

T«T C 3 (O, (9) 

and the quantity 

Tr « (O.ZS^ » (O.ZS^ (f^)" 3 = ffi/10»K)(£|||)) 1/3 (10, 

We take T 9 aver in this form for simplicity in calculating the average core density p. The 
neutrino luminosity can then be expressed through integrating the Lane-Emden function, 

L v w (4 x lO^ergcm-^" 1 ) e 1 e 2 (T C /10 9 K) 9 (47ra 3 ) J €l ^ 2 ^(0^ 

= (4x 10 15 ergcm- 3 s^) ei 6 2 (T^) 9 (-£^-) 3 (4vra 3 ) J* 0§(0d£ 

JO 4 y 3 d « 

= (1.6 x 10 18 erg cm" 3 s- 1 ) ei e 2 (T 9 aver ) 9 (4^/3) (11) 

Apparently most of the neutrinos are emitted near the black hole formation 
point, where Tg ver is the highest. Therefore, to a good approximation we can take 
the radius of the core R to be (3 r s « 3 x 10 10 f3 M^ c cm where (3^1, with a 
characteristic dynamic time of (3 t s « (3 M^ c sec. The volume averaged core density is 
p « M HC /(4vr/? 3 r s 3 /3) = 1.83 x 10 6 /3~ 3 (M 5 HC )" 2 gem" 3 . On the other hand, the entropy 
per baryon in the homologous core (FWW) is 

S?° = a = 2.9 a(M^) 1/2 = 2.9 a(M 5 HC ) 1/2 (M'f 7M 5 HC ) 1/2 (12) 

where a = 'S'ioo/^iooj ^he factor by which the initial entropy is reduced by neutrino emission 
in the course of core collapse. Therefore we conclude that 




From eq. flU) we find a 2 (M* nit /M 5 HC ) « ^5.5/2 « 1.66. Therefore T 9 aver « 13 ((3 2 M^ ) -1 / 2 . 
As a result, eq. flllD becomes 

L„« 2 x 10 60 e 1 e 2 /3- 6 (M 5 HC )- 1 - 5 ergcm- 3 s- 1 (14) 
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An accurate estimate of e\ and e 2 requires a 3-D numerical treatment with full general 
relativity. Given all the other uncertainties inherent in the problem of supermassive star 
formation and collapse, it suffices in this paper to give a "ball-park" estimate of these 
factors using the Newtonian picture. By assuming that all neutrinos are emitted in the 
positive radial direction (not a bad one because most neutrinos are emitted from the central 
part of the core, as seen from figure 1), and that the neutrino emission from the core is cut 
off when R = r s , we estimate that t\ pa 1/40. 

Estimating e 2 is trickier. Because most neutrinos are from the central region of the 
core, neutrinos emerging from the core when the core radius is R are mostly emitted 
~ R/c earlier when the core still has a radius of ~ 2R and is less relativistic. We therefore 
apply the core surface redshift factor e 2 ~ 1 — r s /R = 1 — /3" 1 (a factor of \Jl — r s / R 
from energy redshift and an additional \Jl — r s / R from time dilation to observers) to all 
neutrinos. Although the Newtonian picture simply doesn't apply to the strong field and 
time-dependent situation, we believe that by applying the above e\ and e 2 to eq. ( [LTD we 
can get a "ball-park" estimate for the neutrino luminosity 

L„~5x 10 58 /r 6 (l-/r 1 )(M 5 HC )- 1 - 5 ergcnr 3 s- 1 . (15) 

The time evolution of the neutrino luminosity to an observer at infinity is shown in 
figure 2, with (5\/l — being crudely the observer time axis if the core radius is collapsing 
at the speed of light. The peak luminosity is reached at (3 = R/r s pa 7/6. Integrating over 
time yields the total neutrino energy loss up to the formation of the black hole 

E uloss ~ f°° L u yjl - (3- 1 d((3t s ) ~ 3.6 x 10 57 (M 5 HC )- - 5 ergcm- 3 S - 1 pa 0.036 £ S (M 5 HC )- L5 . 

(16) 

The energy loss through neutrinos cannot be greater than the gravitational binding 
energy itself. Therefore, for M5 ^0.1 the neutrino loss will saturate the limit of the 
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gravitational binding energy and the above scaling with Mf° will break down. In fact when 
^0.1 we do not expect our calculations to apply in the first place because neutrinos 
become trapped in the core instead of freely streaming out. 

So far we have assumed that the radius of the core collapses at the speed of light at 
the black hole formation point. The neutrino energy loss and luminosity can be significant 
larger than we have calculated if rapid rotation or strong magnetic fields slow down the 
collapse significantly. In these cases not only will the time available for neutrino release be 
longer, but the time delay factor e 1 will be less damaging. 

3. The Neutrino Spectrum from the Collapse of a Supermassive Star 

The energy spectra associated with neutrino emission from the annihilation of 
e + e~-plasma with a temperature T is best estimated using a Monte Carlo method: we 
pick an electron-positron pair from a thermal distribution at temperature T, and randomly 
distribute the momentum of one neutrino in all directions in the center-of-mass frame. The 
energies and momenta of the two neutrinos produced are then fixed by momentum and 
energy conservation. We then convert the neutrino energies back into the rest frame, and 
count each annihilation as \M\ 2 events, where the amplitude of the annilihation matrix 
element is (Dicus 1972) 

\M\ 2 = constant • {(a 2 + b 2 ) [(p e - ■ q v ){p e + ■ q v ) + (p e - ■ q p )(p e + • q v )\ + (a 2 - b 2 )m 2 e (q v -q p )}. (17) 

In this equation the p's and g's are four-momenta of e + , e", v and v, and m e is the electron 
rest mass. The coefficients are 

a = l±2sin^, 6 = 0.5, (18) 
where sin 6^ ~ 0.23, and the + sign is for v e , the — sign is for and u T . 
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After repeating the same process many times (several billion times in our case), and 
tallying the total number of events that correspond to various energy bins, we have an 
un-normalized neutrino energy spectrum resulting from e + e~ annihilation in a steady state 
equilibrium plasma. All three neutrino species share a common normalization so that the 
calculation gives a relative flux between z/ e z/ e and v^v^ (or v T v T \ which is around 4.7:1 
as long as T ^ m e m 0.5 MeV. The simulation also yields the temperature dependence 
of the neutrino fluxes, which is proportional to T 8 if T ^ m e , consistent with theoretical 
expectations. 

We fit the neutrino spectrum (normalized energy distribution function) to the form 

fu = T2FM eW-^ + l" (19) 
where the relativistic Fermi integrals are 

POO rph A rp 

F k 0n v ) = / . (20) 

v ' ; Jo e x ~^ + 1 v ; 

We find that for T ^ 0.5 MeV, T v « 1.6T and t] v ~ 2 for all neutrino species. The average 

neutrino energy is ~ 5.5T, higher than that of the ambient e + e~ plasma. This is because 

(1) electrons and positrons with higher energies have larger cross sections for annihilation 

into neutrinos, and (2) the mass of the electron and positron add into the energy of the 

neutrinos. 

Figure 3 shows the arbitrarily-normalized neutrino energy spectrum calculated from 
our Monte Carlo method, and the analytical fit, at T = 1.5 MeV. The analytical fit is 
remarkably good. Similar goodness of fit is also obtained for a variety of temperatures T, 
ranging from 0.5 MeV to 10 MeV, fully covering the typical temperatures inside a collapsing 
supermassive star. 

The neutrino spectrum from the collapse of a supermassive star can be estimated from 
the single-temperature spectrum eq. ( |I9| ) averaged over the entire core and over its time 
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evolution. A precise account is not warranted at this stage. As an estimate, we simply 

take the e + e~ temperature to be that at the peak neutrino emission point in both space 

and time. That is, we estimate it at a position r w 0.14/2 (from figure 1) where neutrinos 

emitted at this position emerge from the core when f3 ~ 7/6 (from figure 2). This turns out 

to be T ~ 0.8 (M^ 10 ) -1 / 2 MeV. Therefore the neutrino spectrum from a supermassive star 

with a core mass of M^ c is to a fair approximation 

1 F 2 
fu « ; 7— , (21) 

(0.8MeV/v/Mp) 3 F ^ 2 ) e (E^F/i.2M C v)-2 + 1 



with an average energy 

(E v ) « 4 (M 5 HC ) - 1 / 2 MeV. (22) 

The average energy can be a factor of 2 higher if the collapse is slowed down significantly 
by rotation/magnetic fields such that the travel time difference is much less important. 



4. Neutrino Background from Supermassive Stars and Prospects for Detection 

If supermassive stars were ever ubiquitous in the universe, they would have left a 
significant neutrino background. Assuming that a fraction F of all baryons once resided 
in supermassive stars, and a fraction / of the associated gravitational binding energy was 
released by neutrino emission, then the total background flux from this neutrino background 
is now 

cf> v ~Ffp b c 3 /(E u ), (23) 

where pb is the baryon density today, and (E u ) is the average energy of the neutrinos from 
the collapse of supermassive stars. About 70% of the flux is v e v e , and 15% of the flux is 
Vfj^n or v T v T . Since, from eq. (|16|) and (p2]), /o.i = //0.1 is usually of order 1, and (E u ) ~ 10 
MeV, we have 
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where £\ is Pb divided by the critical density today and h is the present day Hubble 
parameter in units of 100 km/sec/Mpc. 

If F is 10%, a fraction that is certainly allowed by known constraints, then the flux of 
background neutrinos from supermassive stars can be 100 times higher than the neutrino 
background from supernovae (Totani, Sato & Yoshii 1996; Hartmann & Woosley 1997; 
Malaney 1997). This is because metallicity considerations tightly constrain the fraction 
of baryons once in supernova progenitors to be £ 10~ 3 , while the number and energy 
of neutrinos released per baryon may not differ much in the two kinds of collapse event 
(supernovae vs. supermassive stars). The fraction of higher energy neutrinos 10 MeV) 
in the neutrino background from supermassive stars is, however, very likely to be much 
smaller than that of the supernova neutrino background. This is because neutrinos coming 
from supermassive stars suffer a redshift factor (1 + z) if supermassive stars formed and 
collapsed at a redshift of z ^ 1, making the detection of these neutrinos very difficult. 
For example, if z ^ 3, the average energy of the neutrino background from supermassive 
stars will be only <^ 1.4(M| IC ) _0 5 MeV, far below the threshold of the currently-running 
Super Kamiokande experiment (about 7 MeV). Further complicating the detection of 
supermassive star relic neutrinos is the fact that v e in the background will be hopelessly 
buried in the solar neutrino flux, which is ~ 10 9 cm -2 s" 1 at 1 MeV and ~ 10 6 cm" 2 s" 1 
at PS 10 MeV; and the v e flux in the background under the optimal condition F ~ 10% will 
be comparable to or less than the terrestrial V e background from nuclear power stations, 
which ranges from 10 5 to 10 7 cm" 2 s" 1 at the various sites of neutrino detectors (Lagage 
1985). Clearly, an earth-bound detection of the neutrino background from supermassive 
stars will be extremely difficult, if not impossible. 

As worked examples, we calculate the expected event rate from the supermassive 
star neutrino background in several new-generation neutrino detectors, including the 
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currently-running Super Kamiokande. The fiducial masses, thresholds and detection 
channels of these experiments are listed in table 1. Figure 4 shows the expected event 
rates as a function of z, the redshift at which all supermassive stars existed. We adopt 
a neutrino spectrum of the form in eq. ([19|) assuming M^ c = 1, and we assume 100% 
detector efficiencies above the published thresholds. These assumptions are not entirely 
realistic, but they give good order-of-magnitude estimates. The event rates scale linearly 
with the abundance of supermassive stars (represented by F, the fraction of baryons in 
these stars), and /o.i- From figure 4 it can be seen that Super Kamiokande (Totsuka 1997) 
may potentially be able to preclude more than ~ 10% of baryons ever having been in the 
form of supermassive stars at a redshift of 2 C 1 (although such a conjecture may have 
already been ruled out by not seeing them directly!). 



5. Neutrino Burst from Collapse of a Single of Supermassive Star 

It is interesting to calculate the neutrino flux from the collapse of a single supermassive 
star, and see what the requirement would be on detectors to observe such an event. If the 
collapse occured at redshift z, its neutrino fluence now is 

fM HC c 2 , t , 2 ,, ^Hr^/2/lOMeVx /6000Mpcx2 , 

where the luminosity distance to the star. The average neutrino energy now would be 
(E v )/(1 + z)~ 4(M 5 HC )- 1 / 2 (1 + z)- 1 MeV. 

The duration of the neutrino burst, t, is dilated to M|(l + z) sec, or longer if rotation 
and/or magnetic fields prolong the collapse. As long as z ^ 1, then we can conclude that 
di ~ 6000 Mpc. This burst can be converted roughly into numbers of events in detectors 
by scaling from eq. fl24|) and figure 4. The number of events in the detectors are the 
yearly event rates in figure 4 (with F = 1) multiplied by 10 _7 /o.i(^5 IC ) 3 ^ 2 - For example, 
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a neutrino burst from a non-rotating non- magnetized collapsing supermassive star with 
M| = 1 at z — 1 induces 10~ 6 event in Super Kamiokande. This obviously is impossible 
to detect. The v e flux in the burst from the collapse of one supermassive star when the 
v e has a flux and energy comparable to that of the 8 B solar neutrinos will be completely 
swamped by the solar neutrino flux unless z< 1. The V e flux in this example burst, on 
the other hand, is comparable to the terrestrial v e background from nuclear reactors. The 
supermassive star V e component can only stand out from the nuclear reactor neutrinos 
when z <C 1, when their average energy will be higher than the 3 MeV peak energy of 
neutrinos from reactors. Therefore, it seems that the best chance to detect such a neutrino 
burst from a distant supermassive star is in space, in reactions induced by P e , and with 
extremely large detectors. As a simplistic example, to establish a burst signal requires a 
minimum of about 10 u e p — > ne + events to be detected within a duration of ~ 10 seconds. 
On earth, this requires a 3 x 10 10 ton Super Kamiokande-type detector for a supermassive 
star of Mf° — 1 at z — 1. But in space, there may be a significantly lower v e background 
so that the mass of the detector needed will be substantially smaller. One possibility could 
be an AMANDA (Antarctic Muon And Neutrino Detector Array, Lowder et al. 1996) with 
a threshold of <; 5 MeV built on an inactive icy asteroid or satellite of outer planets. 

The frequency of neutrino bursts from supermassive stars (or the frequency of collapse 
of supermassive stars) could be non-negligible. Assuming that they all form and collapse at 
a redshift z, the frequency of these collapse events as observed now is 

2 3 dr p b (l + zfF 
~ 47rr °^ M ' (26) 

where r is the comoving FRW coordinate of these supermassive stars (with earth at the 
origin), a z is the scale factor of the universe at a redshift z (with ao = 1), to is the age of 
the universe, p h « 2 x 10" 29 g/cm~ 2 Vt h h 2 « 5 x 10" 31 g/cm~ 2 (Tytler & Buries 1997) is the 
baryon density of the universe today, and M is the mass of a typical supermassive star. 
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Since dr/dt = c, the speed of light, and r is of order 6000/i 1 Mpc as long as z > 1, this 
frequency is 

O.SFM^ 1 sec" 1 ~ 3 x itfFM^ 1 /day. (27) 

Therefore, even with F ~ 0.004% (i.e., 0.004% of all baryons were incorporated at one time 
into 1O 5 M supermassive stars), these neutrino bursts would be occurring on average about 
once a day, comparable to the occurrence rate of 7-ray bursts. 

6. Summary 

In this paper we have calculated the neutrino luminosity and spectrum from the 
e + e~ — > vv process during the collapse of a supermassive star. We have estimated that 
for a typical supermassive star with a homologous core mass of 1O 5 M , a few percent 
or higher of the total gravitational energy can be carried away by neutrinos, with an 
average energy ~ 4 to 8 MeV, with both the percentage efficiency and the average energy 
depending on the timescale of the collapse. We have further calculated the flux and energy 
spectrum of the neutrino background from a population of supermassive stars. This has 
been done as a function of the redshift of these stars, their abundance, and their efficiency 
in converting gravitational binding energy into neutrino energy. We found that the resulting 
neutrino background could potentially have a higher flux than the supernova neutrino 
background. However, the average energy of this background flux would likely be lower 
because supermassive stars more likely formed at a high redshift, if they had ever formed at 
all. The expected event rates for this background were calculated for several new-generation 
neutrino detectors, including Super Kamiokande. We showed that Super Kamiokande may 
potentially be able to rule out the possibility that more than ~ 10% of baryons could have 
been incorporated in supermassive stars at a redshift z <^ 1. Finally, we calculated the 
flux and energy of the neutrino burst resulting from the collapse of a single supermassive 
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star, and assessed its detectability. These events are extraordinarily difficult to detect. 
The frequency of such burst events, or equivalently the collapse event rate of supermassive 
stars, is also shown to be possibly significant, easily matching or exceeding the frequency of 
occurrence of 7-ray bursts. 

The authors thank Kev Abazajian and Mitesh Patel for helpful discussions. The work 
is supported by NASA grant NAG5-3062 and NSF grant PHY95-03384 at UCSD. 
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Table 1. New-generation neutrino detectors for which rates are calculated. 





Super Kamiokande 




SNO 




HELLAZ 


Detection Channel 


D e p — ► n e + 


v e d - 


-> ppe~ 


(CC) 


1/ e~ — > e~ 




ve~ — > v e~ 


vd - 


-> cpn 


(NC) 




Threshold 


7 MeV (current) 




5 MeV 




0.1 MeV 




5.5 MeV (planned) 










Fiducial Mass 


22.5 kton 




1 kton 




6 ton 
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Figure Caption 



Figure 1. Neutrino production distribution inside the homologous core (un- normalized). 
Solid curve: the distribution when travel time differences are neglected; dashed line: the 
distribution when travel time differences are considered. 



Figure 2. Time evolution of the neutrino luminosity according to eq. (ITH) , for M^ c = 1. 
The neutrino emission stops when the radius of the core R enters the horizon. 

Figure 3. Solid lines: the spectra of the v e v e and v^v^ (or u T u T ) emission due to e + e~ 
annihilation during the collapse of a supermassive star with an assumed core temperature 
T = 1.5 MeV. The dotted line is the analytical fit to the v e v e spectrum. A common 
arbitrary normalization applies to all curves. 



Figure 4. The event rates from the supermassive star neutrino background in several 
new-generation detectors, as a function of the redshift z of supermassive star collapse, for 
M 5 HC = 1. 
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